EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP-20 11-018 
February 28,2011 



Measurement of // 1// production in 
pp collisions at = 7 Te V 



The LHCb Collaboration ^ 



Abstract 

The production of J/^f mesons in proton-proton collisions at ^/s = ITeV is studied with 
the LHCb detector at the LHC. The differential cross-section for prompt J/y production is 
measured as a function of the J/y transverse momentum pj and rapidity y in the fiducial 
region pj- £[0; 14]GeV/c andy £ [2.0 ;4. 5]. The differential cross-section and fraction of 
J/Y from Zj-hadron decays are also measured in the same pj and y ranges. The analysis is 
based on a data sample corresponding to an integrated luminosity of 5.2pb^^. The mea- 
sured cross-sections integrated over the fiducial region are 10.52 ± 0.04 ± 1 .40l[^J ^b for 
prompt J/y production and 1.14 it 0.01 ±0.16 |J.b for J/y from Z7-hadron decays, where 
the first uncertainty is statistical and the second systematic. The prompt J/y production 
cross-section is obtained assuming no J/y polarisation and the third error indicates the 
acceptance uncertainty due to this assumption. 
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1 Introduction 



Understanding J/\\f meson hadroproduction has been a long-term effort both experimentally 
and theoretically. However, despite the considerable progress made in recent years [1], none of 
the existing theoretical models can successfully describe both the transverse momentum (px) 
dependence of the 7/i// cross-section and the 7/i// polarisation measured at the Tevatron. The 
colour- singlet model (CSM) at leading order in [2] underestimates J/x^r production by two 
orders of magnitude [3], and even more at high pi. Including additional processes, such as 
quark and gluon fragmentation [4] leads to a better description of the pj shape at high pi, but 
still fails to reproduce the measured production rates. Computations performed in the frame- 
work of nonrelativistic quantum chromodynamics (NRQCD), where the cc pair can be pro- 
duced in a colour-octet state [5], can explain the shape and the magnitude of the measured J/^f 
cross-section. However, they predict a substantial transverse component for the polarisation 
of J/^f mesons at large pj. This is in disagreement with the CDF J/x^f polarisation measure- 
ment [6], casting doubt on the conclusion that the colour-octet terms dominate production. 
More recent theoretical studies have considered the addition of the gg — > Jjx^cc process to the 
CSM [7, 8], or higher order corrections in gg — )■ J/^fgg [9] and gg — )• J/^i^ggg [10, 11]. 
With these additions, the discrepancy between theoretical predictions and experimental mea- 
surements significantly decreases. However, the agreement is still not perfect, leaving open the 
question of a complete description of 7/i// hadroproduction. The large rate of 7/ 1// production at 
the Large Hadron Collider (LHC) opens the door to new analyses that extend the phase-space 
region explored so far, such as that recently made by the CMS collaboration [12]. In particular, 
the LHCb detector provides the possibility to extend the measurements to the forward rapidity 
region. 

Three sources of 7/i/a production in pp collisions need to be considered when comparing 
experimental observables and theoretical calculations: direct 7/i// production, feed-down 7/1// 
from the decay of other heavier prompt charmonium states like Xci, Xc2 or V^(25), and 7/1// 
from Z?-hadron decay chains. The sum of the first two sources will be called "prompt 7/i//" in 
the following. The third source will be abbreviated as "7/i// from b'\ 

This paper presents the measurement of the differential production cross-section of both 
prompt 7/ 1// and 7/i// from Z? as a function of the 7/i// transverse momentum and rapidity (y) 
with respect to the beam axis in the fiducial region pr E [0; 14] GeV/c and y E [2.0; 4.5]. The 
effect due to the unknown 7/ 1// polarisation is estimated by providing results for the differential 
cross-sections for three extreme polarisation cases. The analysis of a larger data sample is 
needed to measure the 7/i// polarisation over the kinematic range considered. 

2 The LHCb detector, data sample and Monte Carlo simula- 
tion 

The LHCb detector is a forward spectrometer described in detail in Ref. [13]. The detector 
elements are placed along the beam line of the LHC starting with the Vertex Locator (VELO), 
a silicon strip device that surrounds the pp interaction region and is positioned with its sensi- 
tive area 8 mm from the beam during collisions. The VELO provides precise measurements of 
the positions of the primary pp interaction vertices and decay vertices of long-lived hadrons, 
and contributes to the measurement of track momenta. Other detectors used to measure track 
momenta are a large area silicon strip detector located before a 4Tm dipole magnet and a com- 
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bination of silicon strip detectors and straw drift chambers placed after it. Two Ring Imaging 
Cherenkov detectors are used to identify charged hadrons. Further downstream an Electromag- 
netic Calorimeter system (ECAL, Preshower - PRS - and Scintillating Pad Detector - SPD) 
is used for photon detection and electron identification, followed by a Hadron Calorimeter 
(HCAL). The muon detection consists of five muon stations (MUON) equipped with multi- 
wire proportional chambers, with the exception of the centre of the first station, which uses 
triple-GEM detectors. For the data included in this analysis all detector components were fully 
operational and in a stable condition and the main component of the dipole field was pointing 
upwards. 

The LHCb trigger system consists of two levels. The first level (LO), implemented in hard- 
ware, is designed to reduce the LHC bunch crossing frequency of 40 MHz to a maximum of 
1 MHz, at which the complete detector is read out. The ECAL, HCAL and MUON provide the 
capability of first-level hardware triggering. The second level is a software trigger (High Level 
Trigger, HLT) which runs on an event-filter farm and is implemented in two stages. HLTl per- 
forms a partial event reconstruction to confirm the LO trigger decision, and HLT2 performs a 
full event reconstruction to further discriminate signal events. 

The study reported here uses data corresponding to an integrated luminosity of 5.2pb^^ of 
pp collisions produced by the LHC at a centre-of-mass energy of 7TeV in September 2010, 
with at maximum 1.6 MHz collision frequency. The data were collected using two LO trigger 
lines: the single-muon line, which requires one muon candidate with a pj larger than 1 .4GeV/c, 
and the dimuon line, which requires two muon candidates with px larger than 0.56GeV/c and 
0.48GeV/c, respectively. They provide the input candidates for the corresponding HLTl lines: 
the first one confirms the single-muon candidates from LO, and applies a harder pi selection at 
1.8GeV/c; the second line confirms the dimuon candidates and requires their combined mass 
to be greater than 2.5GeV/c^. The HLT2 algorithm selects events having two opposite charged 
muon candidates with an invariant mass greater than 2.9GeV/c^. For a fraction of the data, 
corresponding to an integrated luminosity of 3.0pb^\ the HLTl single muon line was pre- 
scaled by a factor of five. The trigger efficiency is measured independently for the pre-scaled 
data set and for the rest of the sample, and the results subsequently combined. 

To avoid the possibility that a few events with a high occupancy dominate the HLT CPU 
time, a set of global event cuts (GEC) is applied on the hit multiplicities of each sub-detector 
used by the pattern recognition algorithms. These cuts were introduced to cope with conditions 
encountered during the 2010 running period of the LHC, in which the average number of visible 
interactions per bunch crossing was equal to 1.8 for the data used for this analysis, a factor of 
five above the design value, at a time when only one fifth of the event-filter farm was installed. 
The GEC were chosen to reject busy events with a large number of pile-up interactions with 
minimal loss of luminosity. The average number of reconstructed primary vertices in selected 
and triggered events after GEC is equal to 2.1. 

The Monte Carlo samples used for this analysis are based on the Pythia 6.4 generator [14] 
configured with the parameters detailed in Ref. [15]. The EvtGen package [16] was used to 
generate hadron decays, in particular J/^f and Z?-hadrons, and the GEANT4 package [17] for 
the detector simulation. The prompt charmonium production processes activated in Pythia 
are those from the leading-order colour- singlet and colour-octet mechanisms. The Z?-hadron 
production in PYTHIA is based on leading order 2 — )■ 2 QCD processes: qq — )• q'q', qq' — )■ qq' , 
qq — )■ gg, qg — )■ qg, gg — )■ qq and gg — )• gg. QED radiative corrections to the decay J/y ^ l^^l^ 
are generated using the PHOTOS package [18]. 
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3 J/\\f selection 



The analysis selects events in which at least one primary vertex is reconstructed from at least 
five charged tracks seen in the VELO. 7/1// candidates are formed from pairs of opposite sign 
tracks reconstructed in the full tracking system. Each track must have pi above 0.7GeV/c, 
have a good quality of the track fit (;t^^/ndf < 4) and be identified as a muon by ensuring that 
it penetrates the iron of the MUON system. The two muons are required to originate from a 
common vertex, and only candidates with a probability of the vertex fit larger than 0.5% are 
kept. Some charged particles can be reconstructed as more than one track. Duplicate tracks, 
which share too many hits with another track or are too close to another track, are removed. 

7/ ^f from b tend to be produced away from the primary vertex and can be separated from 
prompt 7/i/A, which are produced at the primary vertex, by exploiting the 7/1// pseudo-proper 
time defined as 

Pz 

where z j/y and zpv are the positions along the z-axis (defined along the beam axis, and oriented 
from the VELO to the MUON) of the 7/i// decay vertex and of the primary vertex; is the 
measured 7/i// momentum in the z direction and Mj/^ the nominal 7/i// mass. Given that b~ 
hadrons are not fully reconstructed, the 7/ 1// momentum is used instead of the exact Z?-hadron 
momentum and the variable provides a good estimate of the Z?-hadron decay proper time. For 
events with several primary vertices (68% of the events), the one which is closest to the 7/i// 
vertex in the z direction is selected. 



4 Cross-section determination 

The differential cross-section for 7/i// production in a given (pj,y) bin is defined as 

d'cT ^ N{J/yf^ii+ii-) 

dydpT =Sf X Ctot X =^(7/i// /i+jU") X X A;?T ' 

where A'^(7/i// — )• /i+/i^) is the number of observed 7/i// — )■ ll^ll^ in bin {pT,y), Ctot the 
7/i// detection efficiency including acceptance and trigger efficiency in bin {pT,y), the 
integrated luminosity, SSiJjy^ —f the branching fraction of the 7/i// — )• /i^/i^ decay 

((5.93 ± 0.06) X 10"^ [19]), and Ay = 0.5 and Apx = 1 GeV/c the y and pi bin sizes, respec- 
tively. The transverse momentum is defined as pi = \J p1 + Py and the rapidity is defined as 

1 E I p 

y = -In where (£',p) is the 7/i// four- momentum in the centre-of-mass frame of the 

2 E — p^ 

colliding protons. 

In each bin of pi and y, the fraction of signal 7/i// from all sources, fjj^t, is estimated 
from an extended unbinned maximum likelihood fit to the invariant mass distribution of the 
reconstructed 7/i// candidates in the interval M^^ G [2.95 ; 3.30] GeV/c^, where the signal is 
described by a Crystal Ball function [20] and the combinatorial background by an exponential 
function. The fraction of 7/i// from b is then extracted from a fit to the distribution. 

As an example. Fig. 1 (left) shows the mass distribution together with the fit results for 
one specific bin (3 < < 4 GeV/c, 2.5 <y < 3.0); the fit gives a mass resolution of 12.3 ± 
O.lMeV/c^ and a mean of 3 095.3 ± 0. 1 MeV/c2, where the errors are statistical only. The mass 
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Figure 1 : Dimuon mass distribution (left) and distribution (right), with fit results superimposed, for 
one bin (3 < pj < 4 GeV/c, 2.5 <y < 3.0). On the mass distribution, the solid red line is the total fit 
function, where the signal is described by a Crystal Ball function, and the dashed blue line represents 
the exponential background function. On the distribution, the solid red line is the total fit function 
described in the text, the green dashed line is the prompt // i/A contribution, the single-hatched area is the 
background component and the cross-hatched area is the tail contribution. 



value is close to the known J/^f mass value of 3 096.916 ± 0.01 IMeV/c^ [19], reflecting the 
current status of the mass-scale calibration; the difference between the two values has no effect 
on the results obtained in this analysis. Summing over all bins, a total signal yield of 565000 
events is obtained. 



4.1 Determination of the fraction of J/\\f from b 

The fraction of J/x^f from b, Fb, is determined from the fits to the pseudo-proper time and the 
invariant mass in each bin of pj and y. The signal proper-time distribution is described 
by a delta function at = for the prompt 7/ 1// component, an exponential decay function for 
the J/x^r from b component and a long tail arising from the association of the 7/1// candidate 
with the wrong primary vertex. There are two main reasons for the wrong association: 

1 . Two or more primary vertices are close to each other and a primary vertex is reconstructed 
with tracks belonging to the different vertices, at a position that is different from the true 
primary vertex position. 

2. The primary vertex from which the 7/1// originates is not found because too few tracks 
originating from the vertex are reconstructed, as confirmed by the simulation; the 7/i// 
candidate is then wrongly associated with another primary vertex found in the event. 

In the first case, the positions of the reconstructed and of the true primary vertices are correlated. 
This category of events is distributed around 1^ = for the prompt component, with a width 
larger than the distribution for correctly associated primary vertices. The contribution of 
these events to the distribution is included in the resolution function described below. 

The long tail is predominantly composed of events in the second category. Since the tail dis- 
tribution affects the measurement of the 7/ 1// from b component, a method has been developed 
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to extract its shape from data. The method consists of associating a 7/i// from a given event 
with the primary vertex of the next event in the 7/ 1// sample. This simulates the position of an 
uncorrelated primary vertex with which the 7/ Xjr is associated. The shape of the tail contribution 
to the signal distribution is then obtained from the distribution of 

^next ^ i^J/V - Mj/y, 
' Pz 

where Zpy ' is the position along the z-axis of the primary vertex of the next event. The pri- 
mary vertex reconstruction efficiency is assumed to be equal for prompt 7/i// and 7/1// from b. 
Given the high primary vertex reconstruction efficiency, 99.4%, the uncertainty related to this 
assumption is neglected. 

The function describing the distribution of the signal is therefore 

_ 1^ 

/signalfe;/p,/b,Tb) =/p5(?.)+/b^+ (1 -/b-/p)/%ilfe), (4) 

where /p is the fraction of prompt 7/1// for which the primary vertex is correctly associated, /b 
the fraction of 7/i// from b for which the primary vertex is correctly associated, Tb the Z?-hadron 
pseudo-lifetime and /Jtaii(^z) the probability density function taken as the histogram shape ob- 
tained from the "next event" method and displayed in Fig. 1 (right). The overall fraction of 

fb 

7/i// from b is defined as Fb = — — . This assumes that the fraction of 7/i// from b in the tail 

/p+/b 

events is equal to the fraction measured with the events for which the primary vertex is correctly 
reconstructed. 

The prompt and b components of the signal function are convolved with a double-Gaussian 
resolution function 

o _ ('z-M)2 1 — - 

/resoludonfeM,^l,52,/3) = ^=^e ^ + e (5) 

The widths of the Gaussians are equal to the event-by-event measurement errors cr, multiplied 
by overall scale factors and 52 to take into account possible mis-calibration effects on a. The 
parameter ji is the bias of the measurement and /3 the fraction of the Gaussian with the smaller 
scale factor. For bins with low statistics, a single-Gaussian resolution function is used. 

The background consists of random combinations of muons from semi-leptonic b and c 
decays, which tend to produce positive ty values, as well as of mis-reconstructed tracks from 
decays in flight of kaons and pions which contribute both to positive and negative values. 
The background distribution is parameterised with an empirical function based on the shape of 
the distribution seen in the 7/i// mass sidebands. It is taken as the sum of a delta function 
and five exponential components (three for positive and two for negative t^, the negative and 
positive exponentials with the largest lifetimes having their lifetimes Tl fixed to the same value), 
convolved with the sum of two Gaussian functions of widths Oi and 02 and fractions j8' and 
(1-/30 



/background ( ) 



(l-/l-/2-/3-/4)5(f,) + 0(f,) I /i^+/2^ I + 



Tl T2 



(6) 



where 0(^7) is the step function. All parameters of the background function are determined 
independently in each bin of pj and y, but for bins with low statistics the number of exponential 
components is reduced. The parameters are obtained from a fit to the distribution of the J/ Xjf 
mass sidebands defined as M^^ G [2.95; 3.00] U [3.20 ; 3.25] GeV/c^, and are fixed for the final 
fit. 



The function used to describe the distribution in the final fit is therefore 



f{h'^fp^fb,fj/xff,IJ-,Sl,S2,P,Th) 

fj/\lf 



fv^ih)+fh—- ®/resolutionfe;M,'5l,'52,/3) + (1 -/b-/p)^tailfe) 



+ (7) 



(1 ^ fj/xjf) /background (^z) 



The total fit function is the sum of the products of the mass and fit functions for the signal and 
background. Four bins of pj and y, which contain less than 150 signal 7/l// events as determined 
from the mass fit, are excluded from the analysis. 

As an example, Fig. 1 (right) represents the distribution for one specific bin (3 < px < 
4 GeV/c, 2.5 <y < 3.0) with the fit result superimposed. The RMS of the resolution function 
is 53fs and the fraction of tail events to the number of 7/1// signal is (0.40 ± 0.01)%. As a 
measure of the fit quality, a is calculated for the fit function using a binned event distribution. 
The resulting fit probability for the histogram of Fig. 1 (right) is equal to 87% and similar good 
fits are seen for the other bins. 



4.2 Luminosity 

The luminosity was measured at specific periods during the data taking using both Van der Meer 
scans [21] and a beam-profile method [22]. Two Van der Meer scans were performed in a single 
fill. The analysis of these scans yields consistent results for the absolute luminosity scale with 
a precision of 10%, dominated by the uncertainty in the knowledge of the LHC proton beam 
currents. In the second approach, six separate periods of stable running were chosen, and the 
beam-profiles measured using beam-gas and beam-beam interactions. Using these results, cor- 
recting for crossing angle effects, and knowing the beam currents, the luminosity in each period 
is determined following the analysis procedure described in Ref. [23]. Consistent results are 
found for the absolute luminosity scale in each period, with a precision of 10%, also dominated 
by the beam current uncertainty. These results are in good agreement with those of the Van der 
Meer analysis. The knowledge of the absolute luminosity scale is used to calibrate the number 
of VELO tracks, which is found to be stable throughout the data-taking period and can there- 
fore be used to monitor the instantaneous luminosity of the entire data sample. The integrated 
luminosity of the runs considered in this analysis is determined to be 5.2 ± 0.5 pb^^. 



4.3 Efficiency calculation 

A simulated sample of inclusive, unpolarised 7/i// mesons is used to estimate the total efficiency 
Ctot in each bin of pj and y. The total efficiency is the product of the geometrical acceptance, the 
detection, reconstruction and selection efficiencies, and the trigger efficiency. It is displayed in 
Fig. 2, including both prompt 7/1// and J/y/ from b. The efficiencies are assumed to be equal for 
prompt 7/1// and 7/i// from in a given (pj^y) bin because neither the trigger nor the selection 
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Figure 2: Total //i// efficiency, as a function of pi in bins of y assuming that are produced unpo- 
larised. The efficiency is seen to drop somewhat at the edges of the acceptance. 



makes use of impact parameter or decay length information. This assumption is confirmed with 
studies based on simulation. 

A correction to the efficiency is applied to take into account the effect of the global event 
cuts described in Sec. 2, introduced during data taking to remove high multiplicity events. The 
effect of such cuts on events containing a 7/i// candidate is not well described by the Monte 
Carlo simulation; it is therefore evaluated from data by using an independent trigger, which ac- 
cepts events having at least one track reconstructed in either the VELO or the tracking stations. 
By comparing the number of such triggered signal J/^f candidates before and after GEC, an 
efficiency of (93 ± 2)% is determined from data. 



4.4 Effect of the J/\\f polarisation on the efficiency 

The efficiency is evaluated from a Monte Carlo simulation in which the J/^f is produced unpo- 
larised. However, studies show that non-zero 7/i/A polarisation may lead to very different total 
efficiencies. In this analysis, the efficiency variation is studied in the helicity frame [24]. 
The angular distribution of the /i+ from the J/^f decay is 

°^ 1 +Aecos + Aerf,sin20cos(i)-|-yl(Asin 0cos20, (8) 

dcos0d0 ^ 

where 6 is defined as the angle between the direction of the /i+ momentum in the 7/i/a centre-of- 
mass frame and the direction of the J/ 1// momentum in the centre-of-mass frame of the colliding 
protons, and is the azimuthal angle measured with respect to the production plane formed by 
the momenta of the colliding protons in the J/^f rest frame. When A,^ = and Xq^ = 0, the 
values Xq = +1, —1,0 correspond to fully transverse, fully longitudinal, and no polarisation, 
respectively, which are the three default polarisation scenarios considered in this analysis. 
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The polarisation significantly affects the acceptance and reconstruction efficiencies. The 
relative efficiency change for prompt J/x^f varies between 3% and 30% depending on pi and y, 
when comparing to the unpolarised case. Therefore, the measurement of the differential prompt 
//!//■ cross-section will be given for the three default polarisations and a separate uncertainty due 
to the polarisation will be assigned to the integrated cross-section. 

Three other polarisation configurations were studied, corresponding to {Xq^X^^Xq^) = 
(-1-1,0, — 1), (0, 1/v^, — 1/2) and (0, — 1/v^, — 1/2); these do not produce variations of the 
measured prompt cross-sections larger than those obtained with the default (±1,0,0) scenarios, 
except in some of the bins with 4 <y < 4.5 where the variations are up to 25% larger. 

The Monte Carlo simulation includes polarisation of 7/1// from b as measured at BaBar 
for and 5+ decays [25]. The simulation shows that the polarisation that the J/xjr acquires 
in b decays is largely diluted when using as helicity quantisation axis the 7/i// momentum in 
the laboratory frame instead of the 7/i// momentum in the Z?-hadron rest frame, which is the 
natural polarisation axis. The effect of the 7/1// from b polarisation on the 7/i// acceptance and 
reconstruction efficiencies is less than 0.5%; therefore, no systematic uncertainty is assigned to 
the J/xjf from b cross-section measurement from the unknown 7/1// polarisation. 

5 Systematic uncertainties 

The different contributions to the systematic uncertainties affecting the cross-section measure- 
ment are discussed in the following and summarised in Table 1. 

Due to the finite pj and y resolutions, 7/i// candidates can be assigned to a wrong pj bin 
(inter-bin cross-feed in Table 1). According to Monte Carlo simulations, the average pj resolu- 
tion is 12.7 ±0.2MeV/c and they resolution is (1.4±0.1) x 10~'^. The effect of the j resolution 
is negligible compared to the bin width of = 0.5. The effect of the pj resolution is estimated 
by recomputing the efficiency tables after smearing the pj values with a Gaussian distribution 
of C7 = 20MeV/c. The maximum relative deviation observed is 0.5% and this is the value used 
as systematic uncertainty for the differential cross-section measurement. The effect on the total 
cross-section is much smaller and is ignored. 

The influence of the choice of the fit function used to describe the shape of the dimuon 
mass distribution is estimated by fitting the 7/i// invariant mass distribution with the sum of two 
Crystal Ball functions. The relative difference of 1% in the number of signal events is taken as 
systematic uncertainty. 

A fraction of 7/i// events have a lower mass because of the radiative tail. Based on Monte 
Carlo studies, 2% of the 7/ 1// signal is estimated to be outside the analysis mass window (M^^ < 
2.95GeV/c^) and not counted as signal. The fitted signal yields are therefore corrected by 2%, 
and an uncertainty of 1% is assigned to the cross-section measurements. 

To cross-check and assign a systematic uncertainty to the Monte Carlo determination of 
the muon identification efficiency, the single track muon identification efficiency is measured 
on data using a tag-and-probe method. This method reconstructs 7/i// candidates in which one 
muon is identified by the muon system ("tag") and the other one ("probe") is identified selecting 
a track depositing the energy of minimum-ionising particles in the calorimeters. The absolute 
muon identification efficiency is then evaluated on the probe muon, as a function of the muon 
momentum. The ratio of the muon identification efficiency measured in data to that obtained in 
the Monte Carlo simulation is convolved with the momentum distribution of muons from 7/ 1// 
to compute a correction factor to apply on simulation-based efficiencies. This factor is found to 
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Table 1: Summary of systematic uncertainties. 



Source 


Systematic uncertainty (%) 


Correlated between bins 




Inter-bin cross -feed 


0.5 


Mass fits 


1.0 


Radiative tail 


1.0 


Muon identification 


1.1 


Tracking efficiency 


8.0 


Track 


1.0 


Vertexing 


0.8 


GEC 


2.0 




1 A 

1.0 


Luminosity 


10.0 


Uncorrelated between bins 




Bin size 


A 1 1 C A 

U.l to f J.U 


Trigger 


1.7 to 4.5 


Applied only to J/ Y from b cross- 


■sections, correlated between bins 


GEC efficiency on B events 


2.0 


fits 


3.6 


Applied only to the extrapolation of the bb cross-section 
b hadronisation fractions 2.0 
^{b J/yX) 9.0 



be 1 .024 ± 0.01 1 and is consistent with being constant over the full J/yf transverse momentum 
and rapidity range; the error on the correction factor is used as a systematic uncertainty. The 
residual misalignment between the tracking system and the muon detectors is accounted for in 
this systematic uncertainty. 

Tracking studies have shown that the Monte Carlo simulation reproduces the track-finding 
efficiency in data within 4%. A systematic uncertainty of 4% for each muon is therefore as- 
signed, resulting in a total systematic uncertainty of 8% due to the knowledge of the track 
reconstruction efficiency [26]. The effects of the residual misalignment of the tracking system 
are included in this systematic uncertainty. 

The selection includes a requirement on the track fit quality, which may not be reliably sim- 
ulated. A systematic uncertainty of 0.5% is assigned per track, which is the relative difference 
between the efficiency of this requirement in the simulation and data. 

Similarly, for the cut on the J/xjr vertex probability, a difference of 1.6% is measured 
between the cut efficiency computed in data and simulation. The Monte Carlo efficiency is 
corrected for this difference and a systematic uncertainty of 0.8% (half of the correction) is 
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assigned. 

The unknown 7/ Xjr transverse momentum and rapidity spectra inside the bins affect the ef- 
ficiency values used to extract the cross-section, because an average value of the efficiency is 
computed in each bin. This effect is important close to the edges of the fiducial region. To take 
into account possible efficiency variations inside the bins, each bin is divided into four sub-bins 
(two bins in pj and two bins in y) and the relative deviation between the bin efficiency and the 
average of the efficiencies in the sub-bins is taken as a systematic uncertainty. 

The trigger efficiency can be determined using a trigger-unbiased sample of events that 
would still be triggered if the 7/ i/a candidate were removed. The efficiency obtained with this 
method in each ipr.y) bin is used to check the efficiencies measured in the simulation. The 
systematic uncertainty associated with the trigger efficiency is the difference between the trigger 
efficiency measured in the data and in the simulation. The largest uncertainties are obtained for 
the high rapidity bins. 

The statistical error on the GEC efficiency (2%) is taken as an additional systematic uncer- 
tainty associated with the trigger. This efficiency is extracted from data as explained in Sec. 4.3; 
it is essentially the efficiency of the GEC on prompt 7/i/a. In the simulation, a 2% difference 
is seen between the prompt 7/i/a and the 7/i/a from b efficiency, which is used as an additional 
systematic uncertainty, applied only to the 7/1// from b cross-section measurement. 

Uncertainties related to the fit procedure are taken into account by varying the central 
value of the prompt 7/1// component, ji, which is found to be different from zero. This shift 
could be due to an improper description of the background for events close to = 0. The 
impact of such a shift is studied by fixing ji at two extreme values, /i = — 3fs and /i = 3fs 
and repeating the fit. The relative variation of the number of 7/1// from b, 3.6%, is used as a 
systematic uncertainty. 

The extrapolation to the full polar angle to obtain the bb cross-section uses the average 
branching fraction of inclusive Z^-hadron decays to 7/i// measured at LEP, i.e., ^{b — )■ 7/i//X) = 
(1.16±0.10)% [27]. The underlying assumption is that the Z?-hadron fractions in pp collisions 
at \/s — 1 TeV are identical to those seen in Z ^ bb decays. However, the b hadronisation frac- 
tions may differ at hadronic machines. To estimate the systematic uncertainty due to possibly 
different fractions, the ^{b — )■ 7/i//X) is computed by taking as input for the calculation the 
fractions measured at the Tevatron [28,29] and assuming the partial widths of B^, 5s and A/, 
to J/xjrX to be equal. The relative difference between the estimates of the branching fractions 
based on the fragmentation functions measured at LEP and at the Tevatron, 2%, is taken as 
systematic uncertainty, which only affects the extrapolation of the bb cross- section. 

6 Results 

The measured double-differential cross-sections for prompt 7/i// and 7/i// from b in the various 
{pT,y) bins, after all corrections and assuming no polarisation, are given in Tables 4 and 5, and 
displayed in Figs. 3 and 4. The results for full transverse and full longitudinal polarisation of 
the 7/i// in the helicity frame are given in Tables 6 and 7, and displayed in Fig. 5. 

The integrated cross-section for prompt 7/i// production in the defined fiducial region, sum- 
ming over all bins of the analysis, is 

a (prompt 7/1//, pT < 14 GeV/c, 2.0 < y < 4.5) = 10.52 ± 0.04 ± 1.40l^;5J ^ib, (9) 

where the first uncertainty is statistical and the second systematic. The result is quoted assuming 
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Figure 3 : Differential production cross-section for prompt J/\j/ function of pj in bins of y , assuming 
that prompt J/y are produced unpolarised. The errors are the quadratic sums of the statistical and 
systematic uncertainties. 
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Figure 4: Differential production cross-section for J/\j/ from b as a. function of pj in bins of y. The 
errors are the quadratic sums of the statistical and systematic uncertainties. 



11 



Table 2: Mean pj and RMS for prompt J/y (assumed unpolarised) and J/\j/ from b. The first uncertainty 
is statistical, the second systematic and the third for prompt J/y the uncertainty due to the unknown 
polarisation. 



Pt range Prompt J/xj/ J/xj/ from b 



(GeV/c) yhin (px) (GeV/c) RMSpT(GeV/c) (px) (GeV/c) RMS (GeV/c) 



0- 


14 


2.0 


-2.5 


2.51 ±0.03 ±0.10+°;°^ 


1.80±0.01±0.04+°;°^ 


3.06 ±0.09 ±0.11 


2.22 ±0.02 ±0.04 


0- 


14 


2.5 


-3.0 


2.53 ± 0.01 ±0.06+g:g^ 


1.74±0.01±O.Oll°;5J^ 


3.04 ±0.02 ±0.05 


2.12±0.01±0.01 


0- 


14 


3.0 


-3.5 


2.46 ±0.01 ±0.02+°;°^ 


1.68±0.01±0.0llj;;!]^ 


2.93 ±0.02 ±0.02 


2.03±0.01±0.01 


0- 


13 


3.5 


-4.0 


2.38 ±0.01 ±0.02+°;°^ 


1.61 ±0.01 ±0.01+°;°! 


2.82 ±0.02 ±0.02 


1.92±0.02±0.01 


0- 


11 


4.0 


-4.5 


2.29 ±0.01 ±0.02+°;°^ 


1.50±0.01±0.01+°;°[ 


2.73 ±0.03 ±0.03 


1.77±0.03±0.01 



unpolarised J/yr and the last error indicates the uncertainty related to this assumption. The 
integrated cross-section for the production of J/yf from b in the same fiducial region is 

a (7/v/ from ;?T < 14 GeV/c, 2.0 < 3; < 4.5) = 1. 14 ± 0.01 ±0.16 ^ib, (10) 

where the first uncertainty is statistical and the second systematic. 

The mean and RMS of the pj spectrum in each y bin are displayed in Table 2. The 7/1// 
mesons from Z^-hadron decays have a mean pj and RMS which are approximately 20% larger 
than those of prompt J/xjr mesons. For each J/yr source, the mean pj and RMS are observed to 
decrease with increasing y. 

Table 3 and Fig. 6 show the differential cross-sections ^ integrated over pj, both for unpo- 
larised prompt J/xjr and J/xjr from b. For the two production sources, the cross-sections decrease 
significantly between the central and forward regions of the LHCb acceptance. 
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Figure 5: Differential production cross-section for prompt as a function of pj in bins of y, assuming 
full transverse (left) or full longitudinal (right) J/y polarisation. The errors are the quadratic sums of the 
statistical and systematic uncertainties. 
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Table 3: ^ in nb for prompt J/\^ (assumed unpolarised) and 7/va from b, integrated over /jt- The first 
uncertainty is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. 



Pi range (GeV/c) 


y 


bin 


Prompt y/i// 


7/V/ from b 


0- 


14 


2.0 


-2.5 


5504 d= 83 ±381 ±726 


697±27±40±96 


0- 


14 


2.5 


-3.0 


5096±21±142±672 


608± 7±13±84 


0- 


14 


3.0 


-3.5 


4460 ± 14 ± 59 ±589 


479± 5± 5±66 


0- 


13 


3.5 


-4.0 


3508 ± 12 ± 40±463 


307 ± 4± 3 ±42 


0- 


11 


4.0 


-4.5 


2462 ± 12 ± 48 ±325 


180± 4± 3±25 




y y 

Figure 6: Differential production cross-section as a function of y integrated over pt;, for unpolarised 
prompt J/^f {left) and from b {right). The errors are the quadratic sums of the statistical and system- 
atic uncertainties. 

6.1 Fraction of //i/a from b 

Table 8 and Fig. 7 give the values of the fraction of 7/i// from b in the different bins assuming 
that the prompt J/^f are produced unpolarised. The third uncertainty in Table 8 gives the devi- 
ation from the central value when the prompt J/ Xjr are fully transversely or fully longitudinally 
polarised in the helicity frame. 

In Fig. 7, only the statistical and systematic uncertainties are displayed, added quadrati- 
cally, but not the uncertainties associated with the prompt 7/i// polarisation. The fraction of 
J/xjf from b increases as a function of pj. For a constant pj, the fraction of 7/i// from b de- 
creases with increasing y, indicating that Z7-hadrons are produced more centrally than prompt 
J/¥- 
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Figure 7: Fraction of J/y from b as a function of pj, in bins of y. 



6.2 Cross-section extrapolation 

Using the LHCb Monte Carlo simulation based on Pythia 6.4 [14] and EvtGen [16], the result 
quoted in Eq. (10) is extrapolated to the full polar angle range 



o{pp bbX) = 



o {J/yr from b, pj < 14 GeV/c, 2.0 < y < 4.5) 
2^{b^J/\j/X) 



(11) 



where a4ji = 5.88 is the ratio of J/y/ from b events in the full range to the number of events in 
the region 2.0 < y < 4.5 and ^{b ^ J/xffX) = (1.16 ±0.10)% is the average branching fraction 
of inclusive Z?-hadron decays to 7/1// measured at LEP [27]. The result is 



g{pp bbX) = 288 ± 4 ± 48 M-b , 



(12) 



where the first uncertainty is statistical and the second systematic. The systematic uncertainty 
includes the uncertainties on the b fractions (2%) and on ^{b J/yfX). No additional uncer- 
tainty has been included for the extrapolation factor a^Ti estimated from the simulation. The 
above result is in excellent agreement with o{pp ^ bbX) = 284 ± 20 ± 49 yih obtained from b 
decays into D^iivX [26]. The extrapolation factor a^Ti has also been estimated using predic- 
tions made in the framework of fixed-order next-to-leading log (FONLL) computations [30], 
and found to be equal to ocj^^^^ = 5.21. 



7 Comparison with theoretical models 

Figure 8 compares the LHCb measurement of the differential prompt 7/ i/a production with sev- 
eral recent theory predictions in the LHCb acceptance region: 
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Figure 8: Comparison of the LHCb results for the differential prompt production for unpolarised 
y/i/A {circles with error bars) with: {top, left) direct Jj\\f production as predicted by LO and NLO 
NRQCD; {top, right) direct J/xjf production as predicted by NLO and NNLO* CSM; {bottom, left) 
prompt J/y production as predicted by NLO NRQCD; {bottom, right) prompt J/y production as pre- 
dicted by NLO CEM. A more detailed description of the models and their references is given in the 
text. 



• top, left: direct J/xjf production as calculated from NRQCD at leading-order in (LO, 
filled orange uncertainty band) [31] and next-to-leading order (NLO), with colour-octet 
long distance matrix elements determined from HERA and Tevatron data (hatched green 
uncertainty band) [32], summing the colour-singlet and colour-octet contributions. 

• top, right: direct production as calculated from a NNLO* colour-singlet model (CSM, 
filled red uncertainty band) [11, 33]. The notation NNLO* denotes an evaluation that 
is not a complete next-to-next leading order computation and that can be affected by 
logarithmic corrections, which are however not easily quantifiable. Direct production as 
calculated from NLO CSM (hatched grey uncertainty band) [7,9] is also represented. 

• bottom, left: prompt J/xjr production as calculated from NRQCD at NLO, including con- 
tributions from Xc and V^(25) decays, summing the colour- singlet and colour-octet con- 
tributions [34]. 
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Figure 9: Comparison of the LHCb results for the differential T/i/A from b production for unpolarised 
y/i/A {circles with error bars) with T/i/A from b production as predicted by FONLL {hatched orange 
uncertainty band). A more detailed description of the model and its references is given in the text. 



• bottom, right: prompt J/\\f production as calculated from a NLO colour-evaporation 
model (CEM), including contributions from Xc and V^(25) decays [35]. 

It should be noted that some of the theoretical models compute the direct J/x^f production, 
whereas the prompt J/x^r measurement includes T/i// from decays and, to a smaller extent, 
y/{2S) decays. However, if one takes into account the feed-down contribution, which has been 
estimated to be of the order of 30% averaging over several experimental measurements at lower 
energies [36], a satisfactory agreement is found with the theoretical predictions. 

Figure 9 shows a comparison of the LHCb measurement of the differential 7/i/a from b 
cross-section with a calculation based on the FONLL formalism [30]. This model predicts the 
Z?-quark production cross-section, and includes the fragmentation of the Z7-quark into Z?-hadrons 
and their decay into //i// mesons. The measurements show a very good agreement with the 
calculation. 



8 Conclusions 

The differential cross-section for 7/ 1// production is measured as a function of the 7/ i/a transverse 
momentum and rapidity in the forward region, 2.0 < y < 4.5. The analysis is based on a data 
sample corresponding to an integrated luminosity of 5.2pb^ collected at the Large Hadron 
Collider at a centre-of-mass energy of y/s = 7TeV, and the contributions of prompt 7/1// and 
7/ 1// from b production are individually measured. The results obtained are in good agreement 
with earlier measurements of the 7/1// production cross-section in pp collisions at the same 
centre-of-mass energy, performed by CMS in a region corresponding to the low rapidity part 
of the LHCb acceptance [12]. This measurement is the first measurement of prompt 7/i// and 
7/ 1// from b production in the forward region at ^/s = 1 TeV. 
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A comparison with recent theoretical models shows good general agreement with the mea- 
sured prompt J/x^r cross-section in the LHCb acceptance at high pi. This confirms the progress 
in the theoretical calculations of J/ Xjr hadroproduction, even if the uncertainties on the predic- 
tions are still large. However, the measurement of the differential cross-section alone is not 
sufficient to be able to discriminate amongst the various models, and studies of other observ- 
ables such as the 7/1// polarisation will be necessary. The measurement of the cross-section for 
y/i// from b is found to agree very well with FONLL predictions. An estimate of the bb cross- 
section in pp collisions at ^/s = 7TeV is also obtained, which is in excellent agreement with 
measurements performed analysing different b decay modes [26]. 
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Table 4: ^j^^ in nb/(GeV/c) for prompt J/xif in bins of the J/\\f transverse momentum and rapidity, 
assuming no polarisation. The first error is statistical, the second is the component of the systematic 
uncertainty that is uncorrelated between bins and the third is the correlated component. 



PT(GeV/c) 


2.Q<y<2.5 


2.5<y<2>.Q 


3.0 < J < 3.5 


0- 


1 


1091 ± 70 ±226 ±144 


844± 13±133±111 


749± 7± 46± 99 


1- 


2 


1495 ± 38 ±282 ±197 


1490 ± 12 ± 39 ±197 


1376 ± 8 ±26 ±182 


2- 


3 


1225 ± 20±109±162 


1214± 9± 24±160 


1053 ± 7± 19 ±139 


3- 


4 


777± 11 ± 44±103 


719± 6± 18± 95 


611± 5± 14± 81 


4- 


5 


424± 6± 22± 56 


392± 3± 12± 52 


325± 3± 9± 43 


5- 


6 


230± 4± 12± 30 


206± 2± 8± 27 


167± 2± 5± 22 


6- 


7 


116± 2± 6± 15 


104± 1± 4± 14 


82± 1± 3± 11 


7- 


8 


64± 1± 3± 8 


57± 1± 3± 7 


44± 1± 1± 6 


8- 


9 


37± 1± 1± 5 


31± 1± 1± 4 


23 ± 1± 1± 3 


9- 


10 


19.3±0.7± 0.5± 2.6 


17.4±0.5± 0.2± 2.3 


12.6±0.4±0.1± 1.7 


10- 


11 


11.6±0.5± 0.3± 1.5 


9.8±0.4± 0.1± 1.3 


7.8±0.3±0.1± 1.0 


11- 


12 


6.7±0.4± 0.2± 0.9 


5.9±0.3± 0.1± 0.8 


4.5±0.3±0.1± 0.6 


12- 


13 


4.6±0.3± 0.2± 0.6 


3.5±0.2± 0.1± 0.5 


2.9±0.2±0.1± 0.4 


13- 


14 


2.9±0.3± 0.1± 0.4 


2.6±0.2± 0.1± 0.3 


1.3±0.2±0.1± 0.2 






3.5<y<A.Q 


4.0<y <4.5 




0- 


1 


614± 6± 23± 81 


447± 5± 28± 59 




1- 


2 


1101± 7± 23±145 


807± 7± 28±107 




2- 


3 


839± 6± 19±111 


588± 6± 22± 78 




3- 


4 


471± 4± 13± 62 


315± 4± 14± 42 




4- 


5 


244 ± 3± 7± 32 


163± 3± 6± 22 




5- 


6 


119± 2± 5± 16 


76± 2± 3± 10 




6- 


7 


59± 1± 2± 8 


34±1.1± 1.4± 4.5 




7- 


8 


29± 1± 1± 4 


17±0.7± 0.8± 2.3 




8- 


9 


15.9±0.5± 0.1± 2.1 


8.5±0.5± 0.4± 1.1 




9- 


10 


8.2±0.4± 0.1± 1.1 


4.1±0.3± 0.2± 0.5 




10- 


11 


4.9±0.3± 0.1± 0.6 


2.2±0.2± 0.1 ± 0.3 




11- 


12 


2.6±0.2± 0.1± 0.3 






12- 


13 


1.2±0.1± 0.1± 0.2 
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Table 5: in nb/(GeV/c) for J/xir from b in bins of the J/xj/ transverse momentum and rapidity. The 
first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. 



PT(GeV/c) 


2.0<y<2.5 


2.5<y<3.0 


3.0<y<3.5 


0- 


1 


107± 23± 22± 15 


75± 4± 12± 10 


60± 2± 4± 8 


1- 


2 


156± 11± 30± 22 


147± 4± 4± 20 


123± 3± 2± 17 


2- 


3 


151± 6± 14± 21 


140± 3± 3± 19 


113± 2± 2± 16 


3- 


4 


105± 4± 6± 15 


98± 2± 2± 14 


75± 2± 2± 10 


4- 


5 


67± 2± 3± 9 


57± 1± 2± 8 


44± 1± 1± 6 


5- 


6 


43± 2± 2± 6 


35 ± 1± 1± 5 


26± 1± 1± 4 


6- 


7 


26± 1± 1± 4 


22± 1± 1± 3 


14.9±0.6±0.5±2.1 


7- 


8 


16.1±0.7±0.8±2.2 


12.1±0.5±0.6±1.7 


9.4±0.4±0.3±1.3 


8- 


9 


10.1±0.6±0.3±1.4 


8.2±0.4±0.8±1.1 


5.3±0.3±0.1±0.7 


9- 


10 


6.5±0.4±0.2±0.9 


5.2±0.3±0.1±0.7 


3.4±0.2±0.1±0.5 


10- 


11 


4.4±0.3±0.1±0.6 


3.2±0.2±0.1±0.4 


2.0±0.2±0.1±0.3 


11- 


12 


3.3±0.3±0.1±0.4 


2.2±0.2±0.1±0.3 


1.5±0.2±0.1±0.2 


12- 


13 


1.9±0.2±0.1±0.3 


1.6±0.2±0.1±0.2 


0.9±0.1±0.1±0.1 


13- 


14 


1.2±0.2±0.1±0.2 


0.9±0.1±0.1±0.1 


0.6±0.1±0.1±0.1 






3.5<y<4.0 


4.0<y<4.5 




0- 


1 


41± 2± 2± 6 


22± 2± 1± 3 




1- 


2 


82± 2± 2± 11 


52± 2± 2± 7 




2- 


3 


71± 2± 2± 10 


42± 2± 2± 6 




3- 


4 


48± 1± 1± 7 


28± 1± 1± 4 




A 

H — 




28± 1± 1± 4 


15.0±1.0±0.6±2.1 




5- 


6 


i< A_i_n7-i-n7-i-0 
iJ.OztU. / ztu. / zt/.z 


y.uzt u. / zt u. J zt 1. J 




6- 


7 


8.6±0.4±0.3±1.2 


5.2±0.5±0.2±0.7 




7- 


8 


5.5±0.3±0.2±0.8 


2.8±0.3±0.1±0.4 




8- 


9 


3.2±0.3±0.1±0.4 


1.5±0.2±0.1±0.2 




9- 


10 


1.8±0.2±0.1±0.2 


0.8±0.2±0.1±0.1 




10- 


11 


1.2±0.2±0.1±0.2 


0.5±0.1±0.1±0.1 




11- 


12 


0.6±0.1±0.1±0.1 






12- 


13 


0.3±0.1±0.1±0.1 
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Table 6: in nb/(GeV/c) for prompt J/xj/ in bins of the J/xj/ transverse momentum and rapidity, 

assuming fully transversely polarised J/y- The first error is statistical, the second is the component of 
the systematic uncertainty that is uncorrected between bins and the third is the correlated component. 



PT(GeV/c) 


2.0<y<2.5 


2.5<y<3.0 


3.0 < J < 3.5 


0- 


1 


1282 ± 83 ±266 ±169 


1058 ± 16±166±140 


924± 9±56±122 


1- 


2 


1751± 44±331±231 


1791 ± 15 ± 47 ±236 


1603± 10± 31±212 


2- 


3 


1438 ± 24±129±190 


1423± 11 ± 28±188 


1182± 7±21±156 


3- 


4 


932 ± 


13 ± 53 ±123 


839± 7± 21±111 


675± 5± 15± 89 


4- 


5 


513± 


7± 27± 68 


455± 4± 14± 60 


358 ± 3± 10± 47 


5- 


6 


278 ± 


4± 15± 37 


238± 3± 9± 32 


184± 2± 6± 24 


6- 


7 


140 ± 


3± 7± 19 


120± 2± 5± 16 


91± 1± 3± 12 


7- 


8 


76 ± 


2± 4± 10 


64± 1± 3± 8 


49± 1± 2± 6 


8- 


9 


44 ± 


1± 1± 6 


34± 1± 1± 5 


25 ± 1± 1± 3 


9- 


10 


23 ± 


1± 1± 3 


19.3±0.6± 0.2± 2.6 


13.7±0.5±0.1± 1.8 


10- 


11 


13.5±0.6± 0.4± 1.8 


10.9±0.4± 0.1± 1.4 


8.5±0.4±0.1± 1.1 


11- 


12 


7.7±0.4± 0.3± 1.0 


6.4±0.3± 0.1± 0.8 


4.9±0.3±0.1± 0.6 


12- 


13 


5.2±0.3± 0.2± 0.7 


3.8±0.3± 0.1± 0.5 


3.1±0.2±0.1± 0.4 


13- 


14 


3.3±0.3± 0.1± 0.4 


2.8±0.2± 0.1± 0.4 


1.4±0.2±0.1± 0.2 






3.5<y<4.0 


4.0<y <4.5 




0- 


1 


728 ± 


7± 27± 96 


530± 6± 33± 70 




1- 


2 


1246 ± 


8± 26±164 


902± 7± 31±119 




2- 


3 


913± 


6± 21 ±120 


631 ± 6± 24± 83 




3- 


4 


505 ± 


4± 14± 67 


334± 4± 15± 44 




4- 


5 


262 ± 


3± 8± 35 


172± 3± 7± 23 




5- 


6 


128 ± 


2± 5± 17 


79± 2± 3± 11 




6- 


7 


63 ± 


1± 2± 8 


36± 1± 2± 5 




7- 


8 


32± 


1± 1± 4 


18.3±0.7± 0.8± 2.4 




8- 


9 


17.1±0.6± 0.2± 2.3 


8.9±0.5± 0.4± 1.2 




9- 


10 


8.8±0.4± 0.1± 1.2 


4.3±0.3± 0.2± 0.5 




10- 


11 


5.2±0.3± 0.1± 0.7 


2.4±0.2± 0.1 ± 0.3 




11- 


12 


2.8±0.2± 0.1± 0.4 






12- 


13 


1.3±0.1± 0.1± 0.2 
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Table 7: in nb/(GeV/c) for prompt J/xj/ in bins of the J/xj/ transverse momentum and rapidity, 

assuming fully longitudinally polarised J/\j/- The first error is statistical, the second is the component of 
the systematic uncertainty that is uncorrected between bins and the third is the correlated component. 



PT(GeV/c) 2.0<y<2.5 2.5<y<3.0 3.0<y<3.5 



0- 


1 


839± 54±174±111 


601 ± 


9± 94± 79 


543± 5± 33± 72 


1- 


2 


1157± 29 ±219 ±153 


1114± 


9± 29±147 


1073 ± 7 ±21 ±142 


2- 


3 


945 ± 16 ± 


84 ±125 


938 ± 


7± 19±124 


865± 5± 16±114 


3- 


4 


583 ± 


8± 


33 ± 77 


559 ± 


4± 14± 74 


514± 4± 11± 68 


4- 


5 


315± 


4± 


16 ± 42 


307 ± 


3± 9± 41 


274± 2± 8± 36 


5- 


6 


171± 


3± 


9± 23 


163 ± 


2± 6± 22 


140± 2± 4± 19 


6- 


7 


87 ± 


2± 


5± 12 


83 ± 


1± 3± 11 


70± 1± 3± 9 


7- 


8 


48 ± 


1± 


2± 6 


46 ± 


1± 2± 6 


38± 1± 1± 5 


8- 


9 


29 ± 


1± 


1± 4 


25 ± 


1± 1± 3 


19.8±0.5±0.1± 2.6 


9- 


10 


14.9±0.5± 0.4± 2.0 


14.5±0.4±0.2± 1.9 


10.8±0.4±0.1± 1.4 


10- 


11 


9.1±0.4± 0.3± 1.2 


8.3±0.3±0.1± 1.1 


6.7±0.3±0.1± 0.9 


11- 


12 


5.3±0.3± 0.2± 0.7 


5.0±0.3±0.1± 0.7 


4.0±0.2±0.1± 0.5 


12- 


13 


3.7±0.2± 0.1± 0.5 


3.0±0.2±0.1± 0.4 


2.5±0.2±0.1± 0.4 


13- 


14 


2.3±0.2± 0.1± 0.3 


2.3±0.2±0.1± 0.3 


1.2±0.1±0.1± 0.2 






3.5<y<4.0 


4.0<y <4.5 




0- 


1 


468 ± 


4± 


21 ± 62 


341 ± 


4± 21 ± 45 




1- 


2 


892 ± 


5± 


18±118 


667 ± 


6± 23 ± 88 




2- 


3 


721 ± 


5± 


16 ± 95 


517± 


5± 20± 68 




3- 


4 


415± 


3± 


12± 55 


282 ± 


4± 13± 37 




A— 




215± 


2± 


7± 28 


148 ± 


2± 6± 20 




5- 


6 


104 ± 


1± 


4± 14 


69 ± 


2± 3± 9 




6- 


7 


51± 


1± 


2± 7 


31± 


1± 1± 4 




7- 


8 


26 ± 


1± 


1± 3 


15.8±0.6±0.7± 2.1 




8- 


9 


13.9±0.5± 0.1± 1.8 


7.6±0.4±0.3± 1.0 




9- 


10 


7.1±0.3± 0.1± 0.9 


3.6±0.3±0.2± 0.5 




10- 


11 


4.3±0.2± 0.1± 0.6 


2.0±0.2±0.1± 0.3 




11- 


12 


2.3±0.2± 0.1± 0.3 








12- 


13 


1.0±0.1± 0.1± 0.1 
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Table 8: Fraction of J/\\f from b (m%) ia bins of the J/\\f ti-ansverse momentum and rapidity. The 

first uncertainty is statistical, the second systematic (uncorrelated between bins) and the third is the 
uncertainty due to the unknown polarisation of the prompt J/\\r; the central values are for unpolarised 



Pt{GcY/c) 


2.0<y<2.5 


2.5 <y<3.0 


3.0<y<3.5 




0- 


1 


8.9 ±1.7 ±0.3+^;^ 


8.2±0.4±0.3+^J 


7.4±0.3±0.3: 


n.3 

-2.5 


1- 


2 


9.4±0.7±0.3+^;^ 


9.0±0.2±0.3t^;^ 


8.2±0.2±0.3: 


n.i 

-2.1 


2- 


3 


11.0±0.5±0.4li;^ 


10.3±0.2±0.4;i;4 


9.7±0.2±0.3: 


M.O 
-1.9 


3- 


4 


11.9±0.4±0.4l^j 


12.0±0.2±0.4;^| 


11.0±0.2±0.4: 


F0.9 
-1.8 


4- 


5 


13.6±0.4±0.5+|^ 


12.7±0.3±0.5;^:^ 


11.9±0.3±0.4: 


M.O 
-1.9 


5- 


6 


15.7±0.5±0.6+^-4 


14.6±0.4±0.5;i;^ 


13.6±0.4±0.5: 


M.l 
-2.1 


6- 


7 


18.4±0.7±0.7+^| 


17.5±0.5±0.6;i;^ 


15.4±0.5±0.6: 


M.2 

-2.3 


7- 


8 


20.1 ±0.8 ±0.7+^-6 


17.6±0.7±0.6;^:^ 


17.8±0.7±0.6: 


M.3 
-2.5 


8- 


9 


21.4±1.0±0.8+^;^ 


21.2±0.9±0.8+^;^ 


18.6±1.0±0.7: 


M.4 

-2.6 


9- 


10 


25.3±1.4±0.9+|^ 


23.1±1.2±0.8+^j 


21.5±1.3±0.8: 


M.3 
-2.5 


10- 


11 


27.6 ±1.7 ±1.0+11 


24.6±1.5±0.9+^j 


20.2 ±1.7 ±0.7: 


M.3 
-2.5 


11- 


12 


32.8 ±2.2 ±1.2+1^ 


27.0 ± 2.0 ± 1.0+^ j 


24.7 ±2.2 ±0.9+^;^ 


12- 


13 


28.9 ±2.6 ±1.0+^;^ 


31.3±2.6±1.1+^J 


24.1 ±2.8 ±0.9+^;^ 


13- 


14 


29.8±3.6± 1.1+^1 


26.5 ±2.9 ±1.0+^ J 


32.5 ±4.1 ±1.2+1 J 



0- 


1 


1- 


2 


2- 


3 


3- 


4 


4- 


5 


5- 


6 


6- 


7 


7- 


8 


8- 


9 


9- 


10 


10- 


11 


11- 


12 


12- 


13 



3.5 <y<4.0 



4.0<y<4.5 



6.3 ±0.3 ±0.2+°;^ 
6.9±0.2±0.2+°J 
7.9±0.2±0.3+?;^ 
9.3±0.3±0.3+?;^ 
10.2±0.3±0.4+°;f 
11.6±0.5±0.4+°j 
12.7±0.6±0.5+?;^ 
15.7±0.9±0.6+};^ 
16.6 ±1.2 ±0.6+10 
18.0±1.6±0.6+^;} 
19.8±2.2±0.7+^;} 
19.5±2.9±0.8+^;^ 
21.9±4.4±0.8+^;^ 



4.8±0.4±0.2;°;^ 
6.1 ±0.2 ±0.2+°;^ 
6.7±0.3±0.2;°;^ 
8.1±0.4±0.3+g;^ 
8.4±0.5±0.3+g;^ 
10.7 ±0.7 ±0.4+°;^ 
13.3±1.1±0.5+?:^ 
13.7±1.4±0.5+?;f 
15.2±2.0±0.5+?;^ 
17.0±2.9±0.6t?;7 
17.7 ± 3.9 ± 0.6+° j 
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